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Abstract. 

The one-dimensional (Id) Anderson model (AM), i.e. a tight-binding chain with 
random uncorrelated on-site energies, has statistical anomalies at any rational point 
/ = where a is the lattice constant and A^; is the de Broglie wavelength. We develop 
a regular approach to anomalous statistics of normalized eigenf unctions ip{r) at such 
commensurability points. The approach is based on an exact integral transfer-matrix 
equation for a generating function <I>r('u, 0) {u and (j) have a meaning of the squared 
amplitude and phase of eigenfunctions, r is the position of the observation point). This 
generating function can be used to compute local statistics of eigenfunctions of Id AM 
at any disorder and to address the problem of higher-order anomalies at / = ^ with 
q > 2. The descender of the generating function Vr{<j>) = ^r{u = 0, 0) is shown to be 
the distribution function of phase which determines the Lyapunov exponent and the 
local density of states. 

In the leading order in the small disorder we have derived a second-order partial 
differential equation for the r- independent ("zero- mode") component ^{u,<f>) at the 
E = {f = anomaly. This equation is nonseparable in variables u and (j>. Yet, we 
show that due to a hidden symmetry, it is integrable and we construct an exact solution 
for ^(u, (f>) explicitly in quadratures. Using this solution we have computed moments 
Im = iV(|-0p™) (m > 1) for a chain of the length N ^ oo and found an essential 
difference between their m-behavior in the center-of-band anomaly and for energies 
outside this anomaly. Outside the anomaly the "extrinsic" localization length defined 
from the Lyapunov exponent coincides with that defined from the inverse participation 
ratio ("intrinsic" localization length). This is not the case at the E = anomaly where 
the extrinsic localization length is smaller than the intrinsic one. At i? = one also 
observes an anomalous enhancement of large moments compatible with existence of 
yet another, much smaller characteristic length scale. 



PACS numbers: 72.15.Rn, 72.70.-|-m, 72.20.Ht, 73.23.-b 
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1. Introduction 

Concepts and methods of the locahzation theory, which counts its origin from the 
seminal Anderson paper [1], have penetrated almost all branches of modern physics[2l|3] 
- from the description of transport in disordered media and the Quantum Hall effect to 
the theory of chaotic systems |1] and turbulence. The one- dimensional (Id) Anderson 
model [1] (AM) - a tight-binding model with a diagonal disorder - is determined by the 
Schrodinger equation for the particle wave function ipi at a site i 

tl^l^i-i + + Sill), = E^pi . (1) 

Here, the nearest neighbor hopping amplitude t is the same for all bonds (below we put 
t = 1); the on-site energy Si is a random variable uncorrelated at different sites and 
characterized by a zero average (^j) = and the variance 

{EiEj) = a^b.j , (2) 

which is a measure of the disorder strength. The variance Eq.([2]) is the only quantity 
that enter our theory at weak disorder <^ 1. For a generic on-site disorder we derive 
equations in terms of the entire distribution function of the on-site energies (uncorrelated 
at different sites): 

^(e) = (5(6-6,)). (3) 

For a finite chain (z = 1, 2, ■ ■ ■ , A^), Eq.([T]) is supplemented with the definition 

^0 = = ?/;jv+i, (4) 

which is equivalent to the hard- wall boundary conditions. In the absence of the disorder 
(e = 0), the normalized wave functions of the chain are = \Jl/N sin(A;j) and the 
corresponding eigenenergies are 

E{k) = 2cos{k), (5) 

where k = tiI/N [l = 1, A^); for an infinite chain (A^ — )■ oo) it fills the interval (0, vr). 
Due to the band symmetry, it is sufficient to consider only k G (0, 7r/2). 

The most studied is the continuous limit of AM, where the lattice constant a — >■ at 
ta^ remaining finite [HI [6l [TJ [H [9l [TOl [11] . There was also a great deal of activity [121 [13] 
aimed at a rigorous mathematical description of Id AM. However, despite considerable 
efforts invested, a lot of issues concerning Id AM still remain unsolved. Among them 
there are effects of commensurability between the de-Broglie wavelength A^; (dependent 
on the energy E) and the lattice constant a, i.e. the anomalous behavior at rational 
values of / = k/n. 

The first known manifestation of commensurability effects was found and described 
quite early [HI |T5] for the simplest objects - the density of states (DoS) and the 
Lyapunov exponent. The latter is defined in terms of a solution of the Schrodinger 
equation ([T]) for a semi-infinite chain {i = 1, 2, ■ ■ ■) supplemented with the definition 
ipQ = aX only one end. For an arbitrary energy E and a generic boundary condition 
ipi 1, the solution to Eq.([l]) is a superposition of two solutions, decreasing and 
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increasing with the increase of i. The increasing part determines the Lyapunov exponent 
7(-E') and the corresponding locahzation length 1{E): 

1 1 1 ^ 

j — 2 

In the continuous model with a weak disorder <^ 1, the Lyapunov exponent ■y^{E) 
and the localization length 

W = (7) 

are smooth functions of energy E. However, according to [HI [15] , for a discrete chain, 
the functions 7(-E) and i{E) possess anomalous deviations from 'j^{E) and io{E) in 
narrow windows of the size oc w around the points k = %/2 (i.e., E = E{k) = 0) and 
k = 7r/3 (i.e., E= E{k) = 1). The Lyapunov exponent sharply decreases at k = 7r/2 
(which is usually associated with an increase of the localization length) but may both 
increase or decrease at k = ir/S depending on the third moment (e?) of the on-site energy 
distribution [1^. It was also conjectured [I5] that progressively weakening anomalies 
(for the weak disorder) may take place at every rational point k/n = m/n (with natural 
m and n) of the band. 

More recently [TTIIIH] it has been found that also the statistics of conductance in Id 
AM are anomalous at the center of band [E = 0, k = tt/2). This is not too surprising, 
because the chain conductance is expressed via the reflection and transition coefficients 
of an electron wave coming from a perfect leads to one of the ends of the disordered 
chain. This problem has a lot of similarities with the problem of calculation of the 
Lyapunov exponent: solutions to both of them are expressed in terms of a probability 
distribution function 7^(0) of the phase 0. The latter can be interpreted as a "phase" 
parameter in the representation of the wave function ipj in the form ipj = aj cos{kj + (j)j) 
with slowly varying real amplitude and phase. The center of band anomaly corresponds 
to an emergent non-triviality of the distribution function V{(j)) as compared to the trivial 
isotropic angular distribution in the continuous problem. 

In a sense, the both problems touch "extrinsic" properties of localization, as the 
Lyapunov exponent describes only the tails of localized wave functions. In the problem 
of the average logarithm of conductance, the extrinsic character of this quantity is set 
by the distance L ^ ^{E) between the ideal leads. In other applications of which 
most important is the interplay between the localization and non-linearity [1], one is 
interested in the number of sites with a high amplitude of the wave function. This 
"intrinsic" picture of localization is better represented by the inverse participation ratio 
(IPR) I J dx (|^/'(x)|^), where ip{x) is a random normalized eigenfunction obeying 
the Shrodinger equation Eq.([l]) with the boundary conditions Eq.(jll) at the both ends 
of the chain. The relation between the two descriptions is not clear yet. In particular, 
it is not known whether the localization length, determined via IPR and sensitive to 
short-range characteristics of eigenfunctions, coincides at points of anomaly k/n = m/n 
with the length determined via Lyapunov exponent Eq.(l6]). This question is one of 
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the motivations of the present paper aimed to develop a formahsm to tackle a class of 
"intrinsic" problems connected with statistics of normalized eigenfunctions. 

Here it is worth noting that IPR and other quantities determined by eigenfuctions, 
are much more difficult to compute than the Lyapunov exponent. The reason is that the 
latter problem does not require a stationary solution to the Shrodinger equation neither 
to obey the boundary conditions at the both ends, nor to be normalized. An attack on 
the "two-end" problem was undertaken in the pioneering paper [T3] but with only limited 
success. The normalization condition imposed on the wave function amplitudes {oj} 
turned out to be very difficult for an analytical treatment. Only few quantities which 
are effectively insensitive to the normalization constraint have been calculated: they are 
the averaged density of states (DOS) and the ratios of wave functions at different sites. 

The difficulty of treating IPR and other quantities determined by local 
eigenfunctions is due to the fact that one needs to deal with an unknown joint probability 
distribution function P(m, 0) of two variables, the phase and the amphtude. 

In this paper we develop a regular approach to treat the anomalies in the 
eigenfunction statistics. The approach is based on a transfer-matrix equation (TME) 
for a generating function of two variables <I>(m,0) {u has a meaning of the squared 
amplitude of wave functions), which is a universal tool to describe properties of a 
generic Id or quasi-ld system. The generating function can be used to compute 
any local statistics of normalized eigenfunctions of Id AM. It also determines a joint 
probability distribution P{u, (f)) of the (squared) amplitude u and phase cf) of the random 
eigenfunction tp ~ i/u cos 0. 

We will concentrate mostly on the study of the principle ( "center-of-band" ) anomaly 
at /c = 7r/2 {E{k) = 0) at weak disorder and show that the corresponding TME for 
$(u,</)) has anomalous terms which make it essentially two-dimensional second-order 
partial differential equation (PDE). This equation is nonseparable in variables u and 
phase 0. Yet, we show that it is integrable and we construct its exact solution explicitly 
in quadratures. 

In the next section [2] we present an elementary derivation of the transfer-matrix 
equation (TME) which is valid at any strength of disorder. This derivation does not 
exploit the supersymmetry method [l9], used in earlier approaches (see [201 El [22] ) . By 
the same token we derive the general expression for the statistical moments of {ipl"^ and 
the mean local density of states in terms of the generating function. In section [3] we 
introduce amplitude-phase variables connected with the representation of eigenfunctions 
in the form ip ~ ^/ucos(f). In these variables, assuming weak disorder cr^ ^ 1, we obtain 
a partial differential equation for the generating function, section HI From that moment 
on we concentrate on the study of the principle, center-of-band anomaly. In section [5] 
we consider a partial differential equation for the case k = 7r/2 {E{k) = 0) and show its 
integrability. Namely, we show that this equation can be factorized in new variables and 
we can construct its solutions. But it turns out that due to non-Hermitian nature of the 
differential operator, there is a continuum of possible solutions. This huge redundancy 
problem is analyzed in section |6] where we show how physical requirements imposed 
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on the generating function allow to find the unique solution. This solution is used in 
section [7] to compute statistical moments of l^/ip. It is shown that these moments cannot 
be described by only one-parameter. This invalidate the one-parameter scaling at the 
anomalous center-of-band point = 0. In Conclusion we point out to the analogies 
with certain dynamical systems and discuss the open problems. 



2. Elementary derivation of the transfer-matrix equation (TME) 

The quantity 

Ur^E) ^^^Y.{\^'r^T'KE-E^)) (8) 

generalizes the concept of the inverse participation ratio (IPR)( given by Eq.Q at 
m = 2) to an arbitrary natural m. Here, is a z/-th eigenfunction of Eqs.([T]) and 
dl]) with the eigenenergy Ey, the summation runs over all states, the angular brackets 
denote the averaging over the ensemble of random site potentials {Sj}, and the quantity 

u{E)=N-^Y.m-E.)) (9) 

V 

is the averaged density of states (DoS) at a given energy E. For a weak disorder cr^ ^ 1, 
the DoS v[E) inside the energy band and outside of "anomalous regions" only slightly 
(by 0(cr^)) differs from the corresponding function Vq{Ei) for the ordered system: 

.(g)«-^o(g)= ^^J|,(^)| . (10) 

However, the difference is known [111 [T5] to be appreciable (~ 0(1)) near the = 
anomaly (see section [3^ . 

For m = 1, the quantity Ylir im=\[Ti E) = N just due to the normalization of wave 
functions. This implies that Im=i{E) = 1 + 0(1/A^), as far from the ends of a long chain 
Im{r,E) is independent of the position r. 

For m > 2, the quantity Im{f, E) can be expressed in terms of the Green's functions 
of the problem, Eqs.([T]) and with the use of a limiting procedure (see, e.g., the review 

my- 

Ur,E) = \im p^{G-f\E^)GrAE-)) ; E^ = E±'-^. (11) 
The Green's functions 

GUE^) = T.^^^^ 



^ E± — Ey 

(with a source at the site r) obey the equations 

{Gj~l^r + G'j+l^r] + ^jGj^r + Sj,r = E±Gj^r (12) 
Co.r = = (jTV+l.r • (13) 

Instead of the supersymmetry approach [19], where Green's functions are represented as 
functional integrals over the usual complex ("bosonic") and Grassmann ( "fermionic" ) 
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variables, here we will present an elementary derivation. The derivation is close in spirit 
to the methods used in refs. [HI [15] . 

For j 7^ r, dividing Eq. (fT2l) by Gj^,- (as is justified below, this quantity differs from 
0), we obtain the recursive equation: 



supplemented with the definitions (see Eq.f lT3|) ) 

l/%^ = = g^, (15) 

Starting with j = 1 and using Eq.f ll4l) to go from j — \ to j, one can find all the q^, 
< j < r as functions of ei, . . . ,er-i- Similarly, starting with j = N and going from 
j to j — 1, one finds all the qf,r<j<Nas functions of £r+i, ■ ■ ■ ,£n- Finally, from 
Eq.f lT2|) at j = r we obtain the quantities of our interest, Gr,r{E±): 



Gr,r{E±) = \t Tj-^ = i dX exp 



±i\ \E± - £r - qr - 



■m 



GZ'' (^+) = r^-^ / dX X"^-' exp 
[m — 2)1 Jo 



i\ \E+ - Er - qr - zr 



Qr-l, 

m>2. (17) 



To justify the transition from Eq. fll2l) to Eq. (fT4l) . we should prove that Gj^E^) ^ 0. 
Note that qi = E + ir]/2 — ei, hence '^qi > 0. Assuming that Qqj^-i > 0, we obtain (for 
j < r): = — $5 (l/g^li) + ri/2 > 0, hence Q^g^ > 0, and by induction this is true 

for any j < r, which means all the corresponding (and Gj,r) are nonzero. Similarly, 
for quantities q~ we find Qqf < 0. And for the case j > r, we obtain qj~ < and 



Qq; > 0. 



3 

As a by-product of this proof, we have confirmed the expected positiveness 
(negativeness) of the imaginary part of the denominator of Gr^r{E+) {Gr^r{E-)), which 
justifies the integral representation in Eq.( fT6l) . 

Using Eq. (ll6l) . the expression Eq. fllll) can be represented in the form (m > 2): 

1 /"OO roo 

I^{r,E) = — lim r/™-i / / dAidA2Ar' e*(Ai-A2)i?-{Ai+A.)V2 

27r(m — 2)! z/(_fc) f?-^+o Jo Jo 

(e-^('^-'^)^'■)7^,_l(Al, A2)^.(Ai, As) . (18) 

Here 

7^,(Al, A2) = (exp [-zAi/g+ + zAa/g/]) ; j < r ; (19) 

^,(Ai, A2) = (exp [-zAig+ + iX2qf]) ; r < j , (20) 

where the averaging in Eqs. fllQI) and fl20|) is performed over random energies £1, . . . , e^-i 
and Er+i, . . . ,eN, respectively. The crucial assumption here is that the random energies 
Er are uncorrelated at different sites. The functions TZj and TZj obey recurrent equations. 
To derive them, we use the following identity for the Bessel function Jo{x): 

/•oo (j 

exp{-z) = - dX'Joi2Vy) — exp{-\'/z) z > 0) , (21) 

JO OA 
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which allows us to convert the factor 1 / qj to qj in the exponent of Eq. 0191] and then to 
apply Eq. ffT^ . As a result, we obtain: 



7^i(Al, A2) = dX[dX', Jo(2yA^) Jo(2yA^) 



[^e'W-A',)(i.^.,)^ e-(^i+^^)''/27^,_l(A;, A',)] (22) 

with the initial condition 7^o(Ai,A2) = 1. Note that the function e~'^(K-K)^j in the 
integrand is statistically independent of quantities ej_i,ej_2, etc. which determine the 
function TZj^i{X[, Ag). In a similar way but proceeding from the site to the site j, we 
derive the recursive equation for TZj^Xi, A2): 



7^i(Al,A2) = dX\dX'^U2^[xJ^f)Jo{2^X2>l,) 



d d r/^,(A;-A',){iJ-.,+,)^e-(^'i+^^)^/2^,+i(A;,A'2)] (23) 



dX[ dX'2 

with the initial condition 'R-NiXi, X2) = 1. For the considered case of site-independent 
statistics of the local disorder, one can see immediately that the function TZN-j{Xi, A2) 
obeys the recursive Eq. fl23l) and equals unity at j = N. Therefore, this function should 
coincide with the function TZj{Xi, A2) and we arrive at an identity: 

7^J(Al,A2) =7^^-J(Al,A2). (24) 

To perform the limit operation rj — )■ +0 in Eq. llTSl) and Eq.( !22l) . we introduce the new 
variables 

s = r/(Ai + A2)/2 ; t; = A1-A2 (25) 

and the new functions: 

Wj{s,v) = njis/r] + v/2,s/r]-v/2) ; TZjis/r] + v/2, s/t] - v/2) = Wn-j{s,v) . (26) 

Using asymptotic form of the Bessel function, integrating by parts in Eq. (l22l) . and 
neglecting infinitely fast oscillating terms, we arrive at the following recursive equation: 



V V 

- + - 

s s' 



e™'^x(^')W^.-i(^',^'), (27) 



where xW) is the characteristic function of the on-site energy distribution J^{e) Eq.(l3]): 
X{v) = J deT{e) = (e"*^'"^ ) . (28) 
In the new variables, Eq. ffTSj) takes the form (m > 2): 



1 



00 roo 



00 



/^(r,E) = ^^^^ J^^,^^^^ / dv i dss^-^e'''''-'xiv)Wr-iis,v)WN-ris,v). (29) 



Free of any limit operation, Eqs. fl27p and are the starting point of our analysis. For 
the Id problem of interest, they were obtained in [22], and still earlier in [20] |2T| in a 
study of the localization transition on the Bethe lattice. The presented here elementary 
derivation of these equations is considerably simpler than the supersymmetry approach 
used in [22l [20] . Also, it allows to establish a relation between the generating function 
Wj{s,v) of our interest and the phase distribution function V{(j)) (see below). 
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3. The "amplitude-phase" variables z and 

3.1. Exact equations in {s,q) and (z,(j)) variables 

To proceed, we introduce the Fourier-transform of Wj{s, v) in the variable v. 

Wj{s,q) = J dve''i''Wj{s,v) . (30) 
The basic equations (I29l).( l271) in the new variables {s,q) take the following form: 

Ur, E) = — -— — / ds s^-^ e-' 

r ^Wr.,{s,q):FiE-q-q')W^.ris,q'), (31) 

2 

Wjis^q) = / dq'TiE - q~' - q') W,.,{sq\q') . (32) 

q J — oo 

Eqs. fl3T]) . fl32|) are exact for uncorrelated on-site energies with the arbitrary distribution 
function J^{e). 

Now we introduce yet another couple of variables, z and (p, determined by: 

s = z cos^ {(p + k) ; q = — txtTT ' (33) 

cos [(p + k) 

and a new function: 

with the boundary condition at j = 

Wo{s,v) = '}Zo{s/r] + v/2,s/r]-v/2) = l 

%{z, <P) = 2vr^(g) 2^jr(^^fc) = - ^/2) • (35) 

Here z G (0, oo); k is determined by the relation E = 2 cos A;; and the "phase" variable 
(f) changes within the interval (0,7r), where there is one-to-one correspondence between 
and q{(f)). Alternatively, we may use so called "extended band" representation, where 
is arbitrary, but the function ^j{z,(f)) obeys the periodicity condition 

$,(^,0 + 7r) = $,(2,0) ; V0, (36) 

and integrations over can be taken over any interval of the length vr. 
In new variables, exact Eq.f l3ip for moments reads 

^mir.E) = — — - / rf0rf0'cos2("'-^)(0) / dzz'^-'e-^'"' 

(m — 2)!z/(i?) JO Jo 

_ / ^ cos cos 0' \ ^ , , , ^ / cos^ ((!) + k) , 

while Eq. (l32|) takes the form: 

$ (^,0) = sin fee- J ^ rci0'^((sinA;(tan0'-tan0)) (^^^,^'-k] . (38) 
cos^ JO \ COS"' y 

From now on the function ^j{z, 0) will be referred to as the generating function which 

determines the statistical moments of distribution. 
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3.2. The Lyapunov and the reflection phase. 



The integrand in Eq.( !37l) for the m-th moment of the quantum-mechanical probabihty 
density Im ~ d'^A^I'") contains cos^™(0)2;™' which suggests the physical meaning of 
and phase and the amplitude of the wave function ip oc y/z cos(0). However, 

one may ask a question what is the meaning of a phase for a wave function in one 
dimensions which may always be chosen real. We answer this question below and show 
that the distribution of phase is related with that of the phase of coefficient r = |r|e*^ 
of reflection from a semi-inflnite disordered chain |25j . 

To this aim, we follow Kappus and Wegner [13] and introduce aj > and (pj 
variables deflned on the link between the sites j (0 < j < r) and j + 1 in such a way 
that (for a flxed site r) 

COS I 

Gj^r{E) = aj cos ; Gj^i^r{E) = aj cos {(f)j + k) =^ q,j = ^— . (39) 

COS 

One can see that such an ansatz is compatible with Eqs. (fT2l) - (fT^ . For the link between 
the sites and 1 we deflne 0o = vr/2, so that 1/go = 0. For brevity, superscripts ± in 
Eq.(l39l) are not indicated. 

Now we use the well-known expression for the Green's function in terms of the 
solutions ■ip'^ and ip'^ to the Schroedinger equation Eq.(Il]) with the arbitrary energy 
E which obey only one of the two boundary conditions: the function ■ip'^ obeying the 
boundary condition ip'^{Q) = 0, while the function tp'^ obeying %p^{N + 1) = 0: 

GrAE)-W^ \V^>(r)#(r'), r > r" ^^^^ 

where We = '^^e^t) V'eI'" + 1) ~ i^Ei^) i'si''^ + 1) is the Wronskian. 
Then Eq. fl39l) is equivalent to 

= v^cos0j ; + l) = ^jCos{(pj + k). (41) 

Eq.f l4ip gives the deflnition of the variables z and (p in terms of the solution to the 
Schroedinger Eq.([T]) with the boundary condition at the left end. As was explained 
above this is exactly the formulation of the problem of Lyapunov exponent. That is 
why we will refer to the phase (p as the "Lyapunov phase". In the absence of disorder 
the relationship Eq.f l4ip is natural, as the shift of phase between the j-th and j + 1-th 
site is indeed equal to k. We will see later on that the so deflned phase has a flat 
distribution V{(p) = ^ far from the end (j ^ ^{E)) of a weakly disordered chain unless 
the energy E is close to the anomalous points k{E) = irp/q. 

The probability distribution V{(p) of the Lyapunov phase is related with the 
distribution Prcf{0) (derived in [25]) of the phase of the coefficient of reflection = 
|r7v|e*^^ from a long chain (of size N ^ oo). Consider a scattering problem deflned by 
Eq.(IT]) on a semi-inflnite (n > 0) chain with no disorder [Sn = 0) at sites n > iV-|- 1 and 
the boundary condition tpn=o = 0. At n > N, the wave function is taken in the form of 
an incident and reflected waves 

= e-^'=('^-^) + r^e*^("-^) ; = e''^^ , (42) 
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the last equation is due to the unitarity. The choice of the "hard wall" boundary 
condition '?/'n=o = is not significant in the limit — oo. The ratios ipn+i/'^Pn = 
with < n < obey the recursive equation 

F„ = ^-£„--^ (43) 

with the boundary condition l/Yn=o = 0. It is crucial that Yn knows only about Ym at 
m < n and is insensitive to the distribution of at m > n. That is why the quantities 
y„ at < < coincide with the corresponding quantities ipn+i/i'n fo'^ semi-infinite 
chain in the Lyapunov problem Eqs. fHOj) and fj^Tl) : 



y„ = ^ = ^^^%i^ (^<n<Ny (44) 

%l)< COS(0„) 

As a consequence, the distribution of random Y^ will be the same as in the bulk of a 
semi- infinite chain and thus governed by the stationary distribution function V{(j)) of 
the Lyapunov phase. On the other hand, it follows from Eq. (H2l) that 

_ ijN+i _ e-^'^ + r^e^ _ cos(^ + k) 

" ^ " l + r, - cos(f ) • ^^^^ 

Comparing Eqs.f l44p and fHS]) one concludes that the distributions of the random 
quantities (pN and On are connected. At A^ — oo the both distributions approach 
their stationary limits, P(0) and Prcf(^), respectively, with the relation: 

PUe) = \rm<t>=e/2. (46) 

Thus the distribution of the Lyapunov phase defined locally at each link of the disordered 
chain and the distribution of the global refiection phase in a semi-infinite disordered 
chain are related in the simplest possible way. 

3.3. Exact relation between the generating function and the phase distribution function 

At vanishing disorder when the on-site energy distribution function J^{e) = S{e), Eq.f l38|) 
equation reduces to 

<l>,+^{z,<p)=e-^'^^'%{z,<p-k). (47) 

In particular, at 2; = using the boundary conditions Eq. (l35|) we obtain: 

<l>,(O,0) = 5(0-jfc-|) . (48) 

As the phase (pj defined in Eq. (14T]) at vanishing disorder is varying like kj with the site 
number j, Eq. (l48|) suggests that ^j{z = 0,0) is the phase distribution function. Now 
we prove that this statement is true at an arbitrary disorder. 
Indeed, from definitions Eqs. ffTOl) and fl26l) we obtain 
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Figure 1. (color online) Stationary distribution of the Lyapunov phase at _E = 
obtained from numerical solution of Eq. (|5ip for Gaussian disorder with the dispersion 
CT^ = 0.1 (red), 1.0 (blue), and 10 (blaek). The dotted line is the analytical solution 
Ea. (|84)) for <t^ which is 7r/2-periodic. With increasing disorder the maxima of the 
distribution move towards (f> ~ ±7r/2 thus breaking the tt/2 periodicity. 



Pi<t>) 




0.0 I 

-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



Figure 2. (color online) Stationary distribution of the Lyapunov phase at fc = 1 
{E = 2cos(l)) obtained from numerical solution of Eq. ([51]) for Gaussian disorder with 
the dispersion ct^ =0.1 (red), 1.0 (blue), and 10 (black). The dotted line is the flat 
distribution V{<j)) = ^ corresponding to cr'^ — > 0. Even at cr^ = 0.1 (localization length 
£q ~ 80) the distribution has a pronounced structure and is far from being a constant 

where it was taken into account the merging of and qf in the hmit 77 —t- 0. Passing 
to the variables z and 0, Eq.f lS^ . we find for the generating function ^j{z, 0), Eq. fl51|) . 
at 2 = 0: 

$,(^ = O,0) = (5(0-0,))=P,(0). (50) 

Thus ^j{z = 0,0) is equal to the phase distribution function Vj{(j)) corresponding to 
Eq. fj^Tl) for solutions to Eq.(IT]) defined on a semi-infinite chain with a generic boundary 
condition 4'j=o = 0, 4'j=i ~ 1- 
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The function Vj{(f)) is of interest in its own right (see e.g. [21 [T5]) as it determines 
the Lyapunov exponent. This function obeys an exact recursive equation 

= rrf0'j^((sinfc(tan0'-tan0)) VM' -k)] (51) 

COS^ JO 

P,=o(0) = 5(0-vr/2), (52) 

which follows immediately from Eq. (l38|) for the generating function $j(2;,0) at 2; = 0. 
Eq.f lSTl) is equivalent to the recursive Eq.(22) in Ref.[Tl]. For weak disorder the 
stationary (site-independent) solution Vj{(j)) — )■ V{4>) can be found analytically ( see 
Ref . |15] and a brief discussion in Sec.5). Examples of the distribution functions 
for different strength of disorder and different energies obtained by straightforward 
numerical solution of Eq. ljSTl) are presented in Fig.l (for E = which corresponds to 
the rational / = | with the smallest denominator) and Fig.2 (for the irrational f = 

The important relation Eq.f lSOj) establishes an exact correspondence between the 
generating function ^j{z = 0,0) and the probability distribution function Vj{(j)) for 
the phase 0j. This identity will be used later for the proper normalization of the 
constructed "stationary" (site-independent) solution $(2;,0). It shows also that the 
generating function $j(2;,0) contains much more information about the system than 
Vj{(j)). This is the reason why the problem of our interest - calculation of the statistical 
moments of normalized eigenfunctions, Eqs.(IH]), ( 1371) . determined by the whole ^j{z, 0)- 
is much more difficult than calculation of the Lyapunov exponent and similar quantities 
determined merely by Vj{(f)). To emphasize the difference between the two problems we 
note that the generalized IPR, Eq. flSTjl . is not linear but bi-linear in ^{z, 0). Hence the 
generating function $(2;, 0) itself cannot be considered as a joint probability distribution 
function of z and 0. The problem of finding the joint probability distribution is not 
simple (see section 1^^ . and involves, as the first step, the calculation of the moments 
Im for integer m > 0. This will be our main goal in this paper. 



3.4- Averaged density of states 



For completeness, we apply the developed formalism to derive an expression for the 
averaged local density of states (LDOS), determined by 



u{E,r) = —Q{GrAE+)). 



TT 



(53) 



Using the representation Eq. (!T6|) for the retarded Green's function, definitions Eqs.( l20|) 
(I20I) . and relation Eq.(l2^. we obtain 
1 



iy{E,r) 



IT 



dX e 



iXE4 



-iXSr 



7^^_l(A,0)7^7v-r(A,0) 



(54) 



Similar to the derivation of Eq. (l371) . we use the definition Eq.( l28ll . relation Eq. (l26l) : 
7^J.(A,0) = Wj{X7]/2,X) = Je-'^''Wj{X7]/2,q)dq/{2TT), integrate over A, and take the 
limit ?7 — > 0, to represent Eq.(j54]) in the form 
dq dq' 



iy{E,r) 



:T{E-q-q')Wr{0,q)WN_r{0,q'). 



(55) 
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In (2;, 0) variables, Eq. (15^ - (IM|) . the above expression reads 

u{E,r) = rd<pd^':F (e - , . - 7''^,. ') $.(O,0)<|.;v-.(O,0') 

Jo \ cos (0 + k) cos (0' + / 

= rd<pd<j>'T(E - , . - "Z"^', ] VM)VM-rm . (56) 

JO \ COS (0 + k) cos (0' + / 

4. Weak disorder. Differential equation for the generating function 

An exact integral equation Eq. fISS]) contains all the information about local statistics 
of eigenfunctions at any uncorrelated on-site disorder. However, only at weak disorder 
the statistical anomalies we are focusing at in this paper are sharp. The point is that 
the region in the energy E (or in the parameter k) where statistics is anomalous is 
proportional to Ai? ~ [15] , and for strong disorder cx^ ~ 1 the anomalies are rounded 
off. That is why in what follows we consider only the case of weak disorder <C 1. 

For the case of a weak disorder, <^ 1, when the "bare" (i.e. the one for the 
continuous model)) localization length £0 ^ 1? the "typical" squared amplitude of 
localized eigenfunctions Ztyp ~ l/i'o << 1 and the exponential factor in front of the 
integral in Eq.f lSSp can be expanded in powers of z. From now we will introduce a 
re-scaled variable 

1 _ 2 sin^ {k) 

-toz 

4 cr^ 

keeping the notation $(n,0) = $(2;, 0)|2=4m/£o ^r the function of this variable. 



u = -ioz = ^z, (57) 



4.I. Expression for moments and DOS in case of weak disorder 

At weak disorder one can replace J^{e) in Eqs. fl37|) and fl56l) by a 5-function and perform 
integration over 0'. Using Eq.([5]) one observes that vanishing of the argument in the 
5-function results in 

COS0' „ , COS0 , ,, ,, , 

777— = 2 COS , , - tan (0' + k) = tan (0 + A;) , (58) 

cos (0' -|- kj cos (0 -|- k) 

from where it follows that 

0' = -0-2A;, mod(7r). (59) 

As a consequence, cos^ {(p + k) = cos^ (0' -|- k), and we arrive at the following expression 
for moments 

^r~l{u,(t) - k)^N~r{u,-(t) - k) . (60) 

Here we have neglected a small (at cr^ ^ 1 ^ £o ^ 1 ) exponent q-^'^^o^'^ i<l>+k)/^o g^j^^j 
used the freedom of shifting the integration interval (as long as it equals to the period 
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27r / d(f)cos^(j)Vr{(p-k)VN^r{-(p-k). (61) 







of the integrand). In a similar way we find from Eq. fl56p the averaged DOS in the weak 
disorder limit: 
iy{E,r) 

ME) 

Up to notations, this relation between DOS and the phase probability distribution Vr{4>) 
coincides with that derived in a different way in [Tl]. For a long chain and far from its 
ends, the function Pr(0) becomes site- independent and the averaged local DOS i'{E, r) 
coincides with ui^E) Eq.(l9]). 

4-2. Joint probability distribution function P{u, (p) 

By definition, the moments Im{r,E) ~ A^d-i/^p™) are expressed via the probability 
distribution function Pr(l'^P) as follows: 

Ur)=N Hj^Priimdm'. (62) 





where the function PrdV^P) is yet to be found. Starting from Eq.( !60l) . we first represent 
it in the form 

l^(^r)=N du c/0[m cos2(0)]™P,(m,0), (63) 
Jo Jo 

which determines a "joint probability distribution function" Pr{u,(j))". This function 
allows to find an average of an arbitrary function fr{u, 0) 

(/,(m,0))=/ du d<Pfr{u,^)Pr{u,<P), (64) 

Jo Jo 

thus providing a more detailed information as compared to P^dV^P)- The connection 
between the two distribution functions is given by: 

^ du rd05(|V^p-McosV)n.(M,0)= r^T-^M^^,0) . (65) 
JO Jo cos^ ycos^fp J 

Now, we will show that the joint probability distribution function Pr{u,(j)) can be 
expressed via the basic object of our study, the generating function $j.(n,0). Using 
in Eq.( l60i) the representation 

m-2 J p+ioo+O Af pt 

(m-2)! du J-ioo+o 2ni t"" ^ ^ 

integrating in u by parts, and changing u — )■ ut, we represent Eq. fl60|) in the form 
Eq.f lB^ . where Pr{u,(j)) is given by 

uo{E) ^ r+ioo+o dt e^*/^o 
Nu{E)u " J-ioc+o 2ni t 
Eqs. (l63il and (l65l) suggest an interpretation of the relation |0^p = ucos^(f) 
as a decomposition of a fast varying (from site to site) eigenfunction ~ 
^Juj cos{(f)j + kj) in terms of two auxiliary slow variables Uj and (f)j. The joint 
distribution function Pr{u,(f)), Eq. (l67|l . allows one also to study separate distribution 
functions of these variables, 

p('^\u)= d(PPr{u,(p) and pW(0) = / duPr{u,(p). (68) 
Jo Jo 



Un(E) r+ioo+0 dt e''*/^° 

Pr{u, 0) = - ')J, du TT- ^ '^-i (^^' <^ - ^) (^^' -(p-k) . (67) 

JSu(E]u J-ioo+o 27n t 
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There is an interesting open question about a relation between the just defined 
distribution function P^'^'^cj)) of the eigenfunction phase for a finite chain with the 
boundary condition = at its both ends, and the distribution function of the 
Lyapunov phase Vr{4>) = ^r{u = 0,0) for a semi-infinite chain. These questions go 
beyond the scope of the current paper. 

What we want to emphasize again at the end of this section is that the complete 
information about the system is conveniently encoded in the generating function 
$j(m,0). Our nearest task is to obtain and solve a differential equation (in the weak 
disorder limit) for this function. 



4-3. Recursive differential operator for the generating function 



For weak disorder [a^ ^1), the function J^{e) in Eq. fl38|) is strongly peaked at e = 
and the integration over 0' is effectively restricted to a narrow vicinity of 0. Expanding 
the remaining part of the integrand in powers of (tan 0' — tan 0) and keeping only the 
first order terms in cr^, we represent Eq. (l38|) in the form: 

4 

+ — mu,(j)) - ci(0) u] 



$j+i(m. 



(69) 



where Ciu, 



is the second order differential operator 

£(m, 0) = C2(0) u'dl + C3(0) {udu - 1) + C4(0) udud^ + C5(0) + c^{^) , (70) 

The coefficients Cj(0) in Eqs. fl69l) and (!70!) are all combinations of cos(20) and sin(20) 
and at first glance do not show any nice structure: 
1, 



[1 + COS(20)], C2(0) = 1 - cos"(20) 



Cl(0) 

C3(0) = -[1 - cos(20) - 2cos^(20)], 
C4(0) = sin(20)[l + cos(20)], 05(0) 



C6(0) 

3 



[1 + cos(20)]^ 



sin(20)[l + cos(20)]. 



(71) 



Note that in the leading order in the disorder strength Eq. (l69|l depends only on the 
variance of the on-site disorder distribution J^{e) through the bare localization length 
£0 given by Eq. fl57|) . 

From Eq. flB^ and the established relation Eq. flSU]) . we can immediately write a 
recursive equation for the phase distribution function 7^(0): 

^ ^' ^^^ V,{<P-k), (72) 



1 + - m 





where 



sin (20) [1 + cos (20)] [1 + cos (20)]^ d 



2 4 9^.- 

Equations (l69l) and ( !72i) are the functional rather than the differential equations because 
of the different phase arguments in the left- and right-hand sides. 
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4-4- Differential equation at rational k/n = m/n ^ 1/2 



In Eqs.( !69l) and ( !72|) the phase argument experiences a finite jump —k at the transition 
from the site j to j + 1. When k = vrm/n with a natural m and n, then after n 
transitions the shift of the phase argument becomes multiple of vr and the functional 
equations can be safely reduced (at weak disorder £o ^ 1) to the differential ones due 
to the periodicity of $(u,</)), Eq.( l36|) . Thus, iterating Eq.f l69|) n times we get a closed 
equation for ^j{u, 0) = ^j{u, cj) — nk)): 



1 + - [£(m,0)-ci(0)m] 

^0 

1 + i [C{u,(t>-{n-l)k) 
to 



Cl( 



[n 



l)k) u] 



^j{u,(l)-nk) . (74) 



Keeping only first order terms in the disorder strength o"^ 



-0 5 



we obtain: 



$j+„(u,0) 



^ 'n—1 m—1 

— C{(j) — r nm/n) — u Ci{(f> — r nm/n) 

-0 Lr=0 r=0 



(75) 



Here the result of summation is extremely sensitive to the particular value of = vrm/n 
and this is the formal reason of an emerging anomaly. Indeed, the functions Eq. (l7T]l 
contain only terms ~ 1, e^^**^, and e^^**^, for which we have 



n-l 



^2i{4>—r nm/n) 



0, 



n-l 



^Ai(ip—r nm/n) 







> 2 
n = 2 . 



(76) 



r=0 r=0 

Thus, for k ^ 7r/2 (i.e., -E'(A;) ^ 0), only ^-independent parts of the coefficients Eq. flTTl) 
survive in Eq. (!75|) . Assuming n <^ Iq, expanding the L.H.S. of Eq. (!75!l . and introducing 
the "continuous" dimensionless coordinate x = 2j/£o along the chain, we obtain: 

3, 



(77) 



The variables u and are separated and one can immediately find a "stationary" (i.e. 
independent of x) solution 

2 



<l>(u,, 



(7J 



This zero mode solution describes the limit of a long chain with the length N ^ io, 
it is the only one which survives at distances x ^ io- This solution has been earlier 
obtained [10] in the continuous limit (n ^ 1). It also arises in the theory of a multi- 
channel disordered wire [191 E]- As follows from Eq.( |77l) . non-zero modes decay at 
distances x = 2j / ~ 1 providing so called "phase randomization": the zero-mode 
solution corresponds to the absolutely isotropic distribution of the phase 0. 

The corresponding moments (?Ti = 1, 2, ...) are found from Eq.( l60|) and are equal 

to: 

jnorm ~ l)' 



(4) 



m—1 



(79) 
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The solution Eq.f l78p corresponds to the following probability distribution of 
squared wave functions {ipl"^ (= zcos^ (</>)) in a long strictly one-dimensional system 
(amazingly, this result was not known before): 

^(l^n^H^ = | '"^y|t''^°^ ^l^l^ (IV^P^o » e-^/^o). (80) 

where N is the chain length. Note that the distribution Eq.l lHOl) is not normalizable, 
as the normalization integral is logarithmically divergent [23]. This divergency is 
an artefact of the zero-mode approximation and is typical to exponentially localized 
wavef unctions. The point is that the eigenfunction statistics changes drastically for 
very small values of the amplitude iV^P^o ^ e~^^^°, where the zero-mode approximation 
no longer applies. This is related with the fact that the envelope (i.e. IV'P averaged 
over oscillations) of the typical localized wave function cannot be significantly smaller 
than lipl"^ ~ e~^^^°. The smaller values of the amplitude \ip\'^ are due to oscillations 
and nodes of wave functions which probability is different from that of the envelope. 
In the exact, iV-dependent distribution function the logarithmic divergency of the 
normalization integral is cut at iV'P^o < e~^^^°- 



5. Center-of-band anomaly, k = 7i/2 



As is seen from Eq. fl76l) for k = irm/n = 7r/2, terms ~ e^^*''^ survive in Eq. flTSj) . This 
leads to a drastic modification of the phase-isotropic equation Eq. fl77l) : 



+ 



L 
3 



u 



$ = 
cos (40) 



[1 - cos (40)] u^df^ + sin (40) udud^ 
dl + 2 cos (40) udu - - sin (40)9<^ - 2 cos (40) 



(81) 



u 



where the differential operator L depends explicitly on 0. The variables u and are not 
separable anymore, which results in an emergent center-of-band {k = -k jl ^ E = 
anomaly: the generating function and the phase distribution function become non- 
isotropic in 0. The variables u and cannot be separated even for the stationary 
variant of Eq.( l8T]) describing the zero mode: 

[L -«]$(«, 0) = 0, (82) 

Yet, due to a hidden symmetry of Eq.l lHT]) . a proper choice of coordinates allows to 
separate variables in the stationary (zero mode) equation Eq.f l82p . This will be done in 
the next subsection. 

For completeness, we conclude this subsection by derivation of an exact expression 
for the stationary distribution function of phase V{(f)) = $(m = 0,0) (see Eq.(l50l)). 
Taking the limit m — > in the stationary variant Eq.( l82|) of Eq. llHT]) . we obtain the 
ordinary differential equation: 







3 + cos (40) ^2 3 . 
4 ^^-2^^" 



3 + cos (40) 



d^- - sin (40) 



- 2 cos (40) 



P(0). 



P(0) 



(83) 
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The only periodic solution to this equation has the form [T5] : 

^^■^"^ i^) = ^ ^ . , (84) 



4^ 



COS (40) 



where the normalization constant provides the equality Jq d(j)V{(f)) = 1. Thus, the 
distribution of the Lyapunov phases in a long weakly disordered chain at the center 
of band {k = 7r/2) is not isotropic but has maxima at (p = According to 

the interpretation Eq. fl39|) of the amplitude-phase variables, this implies the tendency 
towards smaller difference between and I'i/'j+ip, i.e. larger localization length. This 
phenomenon has been coined as "the center- of- band- anomaly " . 

Now we proceed with our much more difficult task: solving not an ordinary but the 
partial differential equation (IHTl) for the generating function $(n,0) of two variables. 



5.1. Hidden symmetry and separation of variables 

The integrability of the stationary equation ( l82l) is shown in three steps. The step one is 
to pass from {u, 0) to a new set of variables (m, v) with v = u cos (20), and to introduce 
a new function 

<l(M,t;) = -$(u,0)|<.os(20)=v« • (85) 
u 

In these variables the zero- mode equation Eq. flSlI) takes a very symmetric form: 

Vu^ - {du Vu^ - f 2 + a„ Vn2 - t;2 a,} $ = ^ $ . (86) 

It is remarkable that the L.H.S. of this equation can be represented as $ where 

the operators Di and belong to the family of three operators from the representation 
of the s/2 algebra: 

Di = a/m2 _ ^2 . D2 = udy + V du ] -D3 = —y/u^ — v"^ (87) 

with the commutation relations: 

[Di, D2] = -Ds, [Ds, Di] = D2, [D2, Ds] = Di. (88) 

Now it is clear that there is a hidden order in a set of coefficients in Eq. (!8T|l resulting 
from the SL(2) symmetry. The latter frequently manifests itself in various scattering 
problems. However, Eq.f lHTl) is connected with even higher symmetry. Introducing a set 
of three additional (mutually commuting) operators 

5i = V, B2 = Vu^ - v^, B3 = u, (89) 

we can represent Eq.f l81l) in the form [Dl + — 3^/2]^ = 0. The operators Di and 
Bi constitute an algebra D Q B with the commutative subalgebra B and commutation 
relations: [D,D] = D (see Eq.(IH8])), [B,B] = 0, and [D,B] = B; in more detail, the 
latter relation looks like: [Di, Bi] = and 

[Di,B2] = Bs ; [D2,Bi]=Bs ; [Ds,Bi] = -B2 ; 

[D,, Bs] = B2 ; [D2, Bs] = B, ; [D3, B2] = B, . (90) 
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It is tempting to interpret D @B as, the algebra of generators of rotations {D) and 
translations [B) of the 3d pseudo-euclidian space R^'^. However, the question of a 
constructive application of this symmetry to the considered problem remains open. 
Below we follow a more prosaic way. 

The next step is to introduce a function 

m{u,v) = {u^ -v'^)-i^{u,v) (91) 

to transform Eq. fl86|) to the Schrodinger-like equation for the function '^{u,v): 

Hm = -{dl + dl) m + U{u, v)m = Q, (92) 
rr/ N 3 + 1;^ 1 u 

U{u,v) = --— ^ + -—2 2- 93 

Finally we introduce the variables 

M + f 2, U — V 2, 

4 = — ^ — = u COS , ?7 = — ^ — = u sin cp . (94) 



It is easy to see that in these variables the operator in Eq. fl92|) splits into two identical 
one- dimensional Hamiltonians 

[H^ + Hr,]^{^,r])=0, (95) 

where is given by: 

Thus, in new variables Eq.f l94p the partial differential equation f l82|) for the generating 
function at = -7r/2 is separable and can be reduced to the two ordinary differential 
equations of the Schrodinger type 

H^MO = A^a(0 ; H.ij^Aiv) = -A^-A(r?) , (97) 

defined on semi-axes ^ > and rj > 0, respectively. The opposite sign of the two 
eigenvalues guarantees the zero-energy solution to Eq. fj95|) . 



5.1.1. Distinctions from the usual quantum mechanics. Although Eqs.f l97j) look like 
ordinary one- dimensional Schrodinger equations on a positive semi-axis, the problem 
we are solving is very different from quantum mechanics. A cornerstone of the latter 
is the Hermiticity of a Hamiltonian which insures that corresponding eigenenergies are 
real and the time evolution of an initial state is unitary. For singular Hamiltonians like 
the one in Eq.f l97p . this property is not given for granted. It requires vanishing the 
boundary term 

coo ^ roo 

Jo 

= mmMo - MOdimMZo = o , (98) 

which arises at integration by parts for any two quadratically integrable functions 

and ip2{0 froKi the Hilbert space. We will show that this condition cannot be fulfilled for 
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the operator Eq. ipBil which eigenfunctions at A 7^ are constructed as a superposition 
of two fundamental solutions with different behaviors at — ?■ 0"*" 

v40~rii + o(0] ; V'+(o~r+[i + o(o]- (99) 

Here the exponents s_ = 1/4 and s+ = 3/4 are the roots of the secular equation 
1) + ^ = 0. For A > the both solutions oscillate at ^ — > oo and are acceptable. For 
A < the both solutions have exponentially decreasing and increasing parts at ,^ — oo, 
therefore only a properly constructed superposition of the two solutions, with cancelation 
of the increasing part, is acceptable. In the special case A = 0, the fundamental solution 
is 

V'olO-e'/'expf-V^V (100) 



while the second solution, ~ ,^^/^exp y—\fijt is unbounded and thus should be ignored. 
Now, considering the equation Eq.( l98|) with the choice 'i^2(0 ~ "^+(0 (f^^ A > 0) 
and ~ ^0(^)5 fi^'i thaX the boundary term vanishes at ,^ — )■ 00 while it is not 

zero at ^ — i- 0: 

This means the Hamiltonian Eq.f l96|) is non-Hermitian in the Hilbert space which 
includes all the (bounded) fundamental solutions. The remedy of the usual quantum 
mechanics [30] to preserve the Hermiticity of singular Hamiltonians like Eq. (p6|l is to 
reduce the Hilbert space to include only one of the two types (characterizing by the 
exponents s±) fundamental solutions. A direct consequence of such a restriction is the 
absence of the bound states of negative energy [30]. Indeed, as mentioned above, for 
A < only a superposition of two fundamental solutions can provide a decrease of the 
wave function at ^ — j- 00. 

There is no such a restriction imposed on the Hilbert space in our problem where 
the eigenvalue A plays an auxiliary role (the zero mode of Eq.f l95p arises as a result of 
cancelation A + (—A) = 0) and does not have meaning of an observable. That is why 
A is allowed to be complex and there is no requirement of Hermiticity of the operator 
Eq.( p6|) . We will see that in contrast to quantum mechanics, bound states (decaying at 
E, — )■ 00) do exist in our problem. Moreover, their spectrum is not discrete but fills a 
sector on the complex A plane. 

5.1.2. Construction of a general solution, A 7^ 0. Solutions iPa{0 with A 7^ can 
be expressed via the Whittaker function W-x^fj,{C,) which obeys the Weber's differential 
equation (see, e.g. Ref. [M], 9.220): 

^ ^-^-(-) + ("i - ^ + ^) ^-U-) = (102) 
and decays at x 00: 
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M_A,^(x) = x^+^e-^/\Fi(^ + ^ + A, 2/i + 1; x) , (104) 

M_A,_^(x) =x^^e-^'/^Fi(^-/i + A,-2/i + l;x), (105) 

where iFi is the confluent hypergeometric function, and 3ft A > 0. For A < 0, the 

equation ( p6|) is mapped on Eq. fll02p by the following identiflcation: 

1 f 1 

A = , ; x = 2v^^ = -V ; /i = -. (106) 

Note that the solution ^^{Cj = W__x,i (^) to Eqs.fl97l). ffT02D at A < 0, which decays 

1 3 

at ^ — )■ oo, contains both a part ~ ^4 and a part ~ ^4 at ^ — 0. This clearly violates 



the condition of Hermiticity Eq.( l98l) and is the reason why the singular Hamiltonian 
Eq.f l96|) does have bound states. Moreover, the spectrum of the Hamiltonian is complex, 
as Eq.f ll06p can be easily extended to complex A with the convention that y/z > at 
z > and has a cut along the semi- axis 2; < 0. From the asymptotic of the Whittaker 
function 

exp(-Fv) i e^oo, (107) 







we flnd the domain on the complex plane of A = |A|e*", where the solution Eq. fll07p is 
decaying: 3ftA > 0, i.e., a G (— 7r/2, 7r/2). A general solution to Eq.(l92l) can be built as 
a superposition of "elementary blocks" 

(108) 

with the following restrictions for Ai and A2: (i) \l + \l = l/(— A) + 1/A = ; (ii) 
3ftAi > 0, 3ftA2 > 0. These inequalities are non-strict because in the case where one of the 
A's is purely imaginary the second one is real due to the condition (i) and this provides 
the required vanishing of the product Eq. (11081) at n — )■ 00. The relation (i) leads to 
the representation Ai = |A|e*", A2 = |A|e*^ with a, /3 G [— 7r/2,7r/2] and /3 = a ± 7r/2. 
Choice of the upper (lower) sign means that Ai is restricted to the IV (I) quadrant of 
the complex plane with a G [— 7r/2,0] (a G [0,7r/2]). 

As a result, the general solution to Eq.( l92i) can be represented as a superposition 



JRA>0,'3A>0 




+ / rf'A c_(A, A) W_,^ i f- W^,,^ i TT ^ V) + ^-(e, V) ,(109) 

J3?A>o,3A<o ■* \4Ay ' 4 V 4A / 

where c±(A,A) are arbitrary functions (restricted by the requirement of convergency of 
the corresponding integrals). 

We mention for completeness, that the Whittaker functions with the second index 
/i = 1/4 which emerges in our problem, constitute a special class. They can be expressed 
in terms of the parabolic cylinder functions Df^{x): i(x) = 2~^{2xY^^D2x-i{V^), 
which reflects the possibility of mapping the problem Eq.f l97p on the problem of two 
"harmonic oscillators" but with an "upside down" potential for one of them. We will 
not explore this correspondence in the present paper. 



Commensurahility effects in Anderson localization 



22 



5.1.3. Solution in the case A = 0. For the eigenvalue A = 0, the mapping Eg. (11061) 
becomes singular. Two eigenf unctions of the operator Eg. (196 p in this simple case 
are elementary functions, increasing and decreasing at ^ — t- oo, respectively. The 
eigenfunction decreasing at ^ — )■ oo is given by Eg. fllOOl) . However, as will be 

shown in the next section, the corresponding solution of Eg. (1971) 



^o(e,r/) = (e^)^/^exp 



;iio) 



does not meet physical reguirements of smoothness of as a function of and 

thus it must be ignored. 

5. i.^. Problem of degeneracy of the general solution. Note that Eg. (11091) possesses 
a huge degeneracy, due to an arbitrary choice of the functions c±(A,A). This is in 
contradiction with an intuitive expectation that the statistics of wave functions in an 
infinite disordered chain should be unigue and independent of the boundary conditions. 
Below we show that the natural physical reguirement of smoothness of as a 

function of helps to determine the solution for the generating function up to a constant 
pre-factor which can be further fixed using the relation Eg. (l50|) and the normalization 
condition Jq V^^'^\(l))d(p = 1 for the anomalous phase distribution function V^"''^\(t)). 

6. Resolving the degeneracy problem and determining the solution for the 
generating function 

6.1. Requirements for the generating function 0) 
The stationary generating function 

<|.(«,0) = {$(e,r/)}^^„_.^^^^„,,„.^=|-|-^vl;(e,r/)| , (111) 

which determines the moments Eg. (l60l) . should obey the following reguirements: 

1. It should vanish at n — t- oo. 

2. It should be periodic in (with the period it). Moreover, for the considered case 
{k = IT /2) the coefficients of Eg.l lHTl) are periodic functions with the period 7r/2 and we 
impose the reguirement $(m, + 7r/2) = $(m, 0) on the stationary solution, too. 

3. It should obey the relation $(u = 0, 0) = V^"'^\(j)), which follows from Egs. (j50l) 
and dHl- 

4. $(n, 0) should be a smooth function of together with all the derivatives with 
respect to 0. It should have no jumps, cusps, etc. 

The first reguirement has been fulfilled due to the proper choice of the integration 
domains in Eg. (ll09p on the complex plane A. The second reguirement is eguivalent 
to \E'(,^,?7) = \E'(?7,^) and can be achieved by the symmetrization of the integrand in 
Eg. (ll09p with respect to the replacement ^ -H- 1]. The third reguirement will be fulfilled 
by the properly chosen behavior of the functions c±(A, A) at |A| — j- (see subsection 16.31 
and Appendix A). Accounting for the fourth, extremely important reguirement is not 
as simple. It will be postponed till subsections 16.31 16.41 
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6.2. From plane to contour integral 

Similarly to the usual coherent states, the set of partial solutions Eq.f llOSp is 
overcomplete, and actually the 2d integration domains in Eq.( ll09p (i.e., the first and 
the fourth quadrants of the complex plane A) can be reduced without losses to a 
Id integration contour. Consider, for instance, the first term in Eq.f llOQp . Note 
that Fj^{\]^,rf) = W^_a i (^) i (Ia) ^ holomorphic function of A in the first 
quadrant, i.e., F+(A;^,?7) depends only on A = pe*" (p = |A|) but not on A = pe~*". In 
polar coordinates (p, a) we have: 

*+(^,^) = I ' da j pdpc+{p,a)F+{pe''';^,r]) = 

Jo JFn 



r^'daf pdpc4p,a)F+{pe"';^,7]) (112) 

Jo JVa 

where the integration contour Tq coincides with the semi-axis A > 0, and the second 
equality is the realization of the possibility to rotate the contour Fq, provided that the 
first argument pe*" of the function remains in the first quadrant. For instance, if 
we choose the contour F^ as a ray which corresponds to the rotation of Fq by the angle 
—a, the variable p on Fq, is represented as |p|e~*". Changing, at a fixed a, the variable 
p in the internal integral: p = te~*", where real t runs from to +oo, we arrive at: 

/•7r/2 r+oo 

^+(^,r/)=/ da dte-^'''tc+{te-''',a)F+{t;^,ri), (113) 
JO Jo 

where we have changed the order of integrations, introduced a new weight function 

C+{t) = t r^\ae-^'"c{te-'",a), (114) 
Jo 

and changed back the notation t -^■ A. Thus without loss of generality we have expressed 
the double integral over the first quadrant in terms of a contour integral 



vl>+(e,r/)= / rfAC+(A)F(A;e,r/). (115) 

Note again, that the key condition for this transformation is a holomorphic dependence 
of F+(A;^,?7) on A. 

Another choice of the integration contour can make the expression more symmetric. 
Namely, rotating the contour Fq by the angle 7r/4, so that A — >■ lAje*"^/^, and introducing 
a new real variable A' by A = A'e*'^/^, we obtain (omitting the prime and re-defining the 
arbitrary function C+(A)) 

t4(„i)-fd\c4X)w_„. [^)w_,,, {^) , (116) 

Here e = e*^/^ e = 6"*^/^ 

Following the same route we obtain a contour integral representation for the second 
term in Eq. ffTOQj) : 

*-(«•") = 1°° 'i>'C-WW.„. (I) IV ... (g) (117) 
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;ii8) 



with C_(A) related with c_(A, A). 

Eqs. flllGp and flll7p determine the generating function $(n,0) Eq. flllip : 

i + v 
{ivY" Jo 

Here we took C+(A) = C„(A) = C(A) to make the integrand symmetric with respect 
to the permutation ^ O 77 in order to fulfill the formulated above requirement 2; the 
function C(A) is a real (without loss of generality) function yet to be determined. Up 
to now we have used only the following loose assumptions on its properties: 

1°. C(A) has no singularities in the first quadrant of the complex A plane; 

2°. at |A| — !■ 00, the integrand in Eq. flllSp decays faster than 1/A; this justifies 
rotations of the contour neglecting contributions of distant arcs. 

6.3. Equation for C{\) 

We begin by determining the behavior of the function C(A) at A — ?■ 0. To this end we 
note that according to the relations Eqs.f lSOp and fl8^ . the generating function $(^,77) 
Eq.f lllSp must tend to a finite limit as ^ — )■ and 77 — 0. 

Re-scaling in Eq. flllSp the integration variable A — mA, we find at m — )■ 



7/3/2 „oo 

'^(O''^)-! ^ ■ ^11/2 / d\C{u\) 

I COS0 sm 01^/^ JO 




1 I I Wn 1 I I + CO. 



(119) 



4 V 4A ; V 4A 

To provide a finite value of the expression (11190 in the limit of vanishing u, we should 
require that (see Appendix A): 

C{X) = , (7(0) = const . (120) 

A2 

A crucial role in further restricting the possible choice of the function (7(A) is played 
by the requirement of smoothness of $(m, (p) as a function of (p (requirement 3 of the 
previous subsection). The generating function Eq.f lllSp is periodic in with the period 
I and it is continuous at the end points = (i.e. 77 = 0), and (p = 7i/2 (i.e. = 0) of 
the interval of the periodicity (0, 7r/2). This is guaranteed by the c.c. term in Eq.f lllSp . 
equivalent to the permutation ^ o 7/. What is not automatically guaranteed is that 
$(^,7/) is smooth as a function of at = and = 7r/2; the smoothness implies the 
continuity of all the derivatives. Amazingly, the requirement of smoothness is sufficient 
to determine the function (7(A) up to a constant pre-factor. As we will see this happens 
because of the special property of the solution Eq.f lllSp encoded in the certain identity 
for the confiuent hypergeometric functions in Eqs. fll03p - fll05p . 

Consider, for instance, the behavior of $(n, 0) at — t- 0, (i.e. 7/ — )■ 0, while ^ ^ u). 
A discontinuity of derivatives at = may arise from the branching of the expression in 
Eq.f lllSp at small rj. Indeed, according to the representation of the Whittaker function 

Eq.f llOSp in terms of M-functions Eqs.f ll04p and f llOSp . we see that 

1 

w^-A,i {^) = {^y[fii\v) + Vvf2{x,v)] , (121) 
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where /i and /2 are analytic functions of rj in the vicinity of 77 = 0. The common factor 
?7^/^ is canceled by the pre- factor in front of the integral in Eq. flllSp . The first term 
in the square brackets of Eq.f ll2ip is regular in the vicinity of 77 = 0, while the second 
one ~ y/r] ~ |0|, is not analytical at 77 = 0. As such a non-analytical behavior is in 
conflict with the requirement 3 of smoothness (section l6.ip . the corresponding part of 
the solution must identically vanish. Extracting this singular (oc ^ in the domain 
V < P^^t ^sing{(,,'n) of general solution Eq.f lllSp . we obtain: 

^ Ua 



C(A) 



W 



M 



er] 
4A 



-Ae- 




0, 



which should be fulfilled for any rj < ^. Vanishing of ^sing{^,v) is equivalent to the 
homogeneous integral equation for the real weight function C{X) with the boundary 
condition C{\ — )■ 0) — > const. 

Similarly, the presence of non- analytical terms and in the special solution 
Eq.f lllOp is the reason why this solution should be ignored. 

The integral equation Eq.f ll22p imposes severe constraints on the function C(A), 
because Eq.f ll22p must be satisfied for arbitrary rj and ^ (at 77 < That is why 
the requirement of smoothness lifts a huge degeneracy and arbitrariness in the possible 
choice of C(A). The existence of even a single (non-zero) solution for C(A) is not evident. 
We will show below that the solution to Eq. fll22p does exist and is unique up to the 
constant pre-factor. 



6.4- Solution for C{\) 

Rotating the integration contours independently for each of the two terms in the 
integrand of Eq.f ll22p and changing A — )■ te and A — > te, respectively, one can make 
the Whittaker function real and take it out of the square brackets. Thus, the integral 
equation Eq.f ll22p takes the form 

'C{et)e^v{-^^) ^ /3 



r( 



it) 



it. 



3 
2 



17] 



C{et)exp{zg) 

r(i + ^t) 



li^i 



it. 



3 
2 



—trj 
IT 



0,(123) 



where the dependence of the integrand on C, and 77 is factorized. The only possibility 
to satisfy this equation for arbitrary ^ and t] is to require the square bracket to vanish 
identically. 

The crucial observation for the possibility to fulfil this condition is an identity for 
the confluent hypergeometric functions [2] 



-2/2 



3 



3 l-Tl 



it. 



(124) 
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With the help of Eq.f ll24p . we find that the square bracket in Eq. (11231) vanishes 
identically for all t] if and only if the function C{t) obeys (for positive t) the condition 

C{et) C{et) 



(125) 



T{\-it) T{\ + it)- 
Now one can immediately guess a solution for C(A): 

C'o(A) = rQ + 6A) rQ + 6A) . (126) 

It is easily seen that the solution Eq. fll26p obeys the both conditions formulated at the 
end of the subsection 16.21 The function Co(A) = C'o(A)/A^/^ (see definition ( I120p ) is an 
analytical function in the domain of our interest (3ft A > 0). Though Co(A) grows at 
|A| — )■ oo, one can check that the integrand in Eq. (11180 decays as l/|Ap for |A| — j- oo. 
This provides the convergence of the integral Eq. flllSI) at large A and justifies rotations 
of integration contours neglecting contributions of infinitely remote arcs. Note also that 
C*o(A) is real at the semi-axis A > 0. 

Now, looking for a general solution to Eq. (11250 in the form 

C'(A) = Co(A)5(A) , (127) 

we obtain the following functional equation for ^(A) at A > 0: 

S{e\) = S{eX). (128) 

The function S'(A) is also supposed to be an analytic function at 3ftA > 0. Equation 
(11281) requires 5'(A) to be an analytical function of z = A^: 

oo 

5(A) = 5:s„A^". (129) 

n=0 

Therefore, being regular in the domain 3ft A > (it is even sufficient to require analyticity 
within a sector |arg(A)| < '/r/4), the function S'(A) must be regular on all the complex 
plane A, i.e. it should be an entire function. Now we apply the condition of convergence 
of the integral over A in Eq. (I118P at large A to find the allowed asymptotic behavior of 
S{X) at |A| — )■ oo. Substituting Eq. (ll27p into Eq. (lll8p and using the asymptotics of the 
Whittaker and F-functions we find that the integrand behaves as A~^S'(A) at A — ?■ oo. 
This means that |5'(A)| should increase not faster than A^. There is only one such entire 
function with the structure of Eq. fll29p : this is a constant iS'(A) = sq = const. Thus 
we have proven the uniqueness of the solution Eq. (ll26p up to a constant factor. This 
factor has to be determined from the relation Eq.(l50l) and the normalization condition 
for the phase distribution function "^(0). Now we may write down the solution for the 
anomalous (at the center of the band) generating function $(""■)(«, 0) in the final form : 



. , roo r (i + eA) r (i + eA) 

214(1) I cos sin 01 1/2 Jo 



A3/2 




(130) 
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Figure 3. (color online) The function 0), in the range u £ [0, 1], </> € [0, 7r/2]. 



where ^ = u cos^ cf) 



, r] = usm' 



2 



e = 




e ^'^Z^. In the appendix we 



demonstrate exphcitly that the obtained solution Eq.f ll3Up does obey the relation 
${an)(^^ = 0,0) = P('^")(0), where p("")(0) is given by Eq.([84]). The 3D plot of the 
function $(m,0) is given in Fig.l. In the next section we apply the solution Eq. fll30p 
for studying moments of the wave function distribution. 

7. Moments of the wave function distribution 

The exact expression Eq.f ll30p for the anomalous (at the center of the band) generating 
function is our main analytic result. It determines statistics of wave function distribution 
at the center-of-band anomaly. Although extensive physical applications go beyond the 
framework of the present work, here we briefly discuss the applicability of the one- 
parameter scaling description for the anomalous statistics. 

As has been mentioned in the introduction, the Lyapunov exponent ■j{E), Eq.([6]), 
sharply decreases in a narrow vicinity of the band center E = [HI [15]. Being dependent 
only on the phase distribution function, the Lyapunov exponent can be easily calculated 
using Eq. ( 184|) at the center-of-band anomaly E = and the trivial homogeneous phase 
distribution 'p("°^™-) (0) = I/tt close to but outside the anomalous region. The ratio of 
real parts of the two corresponding Lyapunov exponents ■janiE = 0) and 'jnormiE ^ 0) 
is given by [HI [15] : 



According to Eq.([H]), this can be interpreted as an increasing localization length at 
the anomaly: 4 ^ ^Tn = 1-094 4. Note that the localization length t""^ = 1/(3^7) 



3f?7,„(E = 0) 

3^ ^norm ( ~ ) 



r[l + cos(40)]p('^")(0) 




0.9139. (131) 
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defined via the Lyapunov exponent 7 characterizes the exponentially decaying tails of 
the localized wave function and for this reason will be referred to as the "extrinsic" 
localization length. In contrast to that we consider the "intrinsic" localization length 
defined via the moments of the inverse participation ratio 

1 



(m-1)! (fnt)m,~i' (132) 

whenever is independent of m in a sufficiently wide interval of m. This localization 
length characterizes the "body" of the localized wave function. 

Comparing Eqs.( ff9l) and ( 1132^ one concludes that away from E = Q anomaly the 
extrinsic and intrinsic localization lengths coincide and are both equal to £0 given by 
Eq.(I7]). To study a relationship between them at the E = Q anomaly, we analyze 
the behavior of moments Im {m > 2) of the anomalous (at the center of band, E = 0, 
k = tt/2) wave function distribution. With the definition Eq.f l57|) . the expression Eq.f l60|) 
in the limit of a long chain takes the form: 

I^-\E = 0) = ^,„,Z!, / d<j> cos'-{<j>) / cinu'"-2[<^>('^")(u,0)]^ (133) 



where the ratio z/o(0)/i^(0) is given by Eq. (l6T]l : 

^ = 47r r^'d<P cos^(0) [$(-)(O,0)]l (134) 
1^0 (0) Jo 



One might expect that the behavior of anomalous moments is similar to Eq.( !79|l 
but with the localization length £q replaced by some other length scale ii^t- This would 
be the scenario of one-parameter scaling which appears to fail at the band center. 

A convenient way to present the results is to plot the reduced moments Rm = 
I^'^^E = 0)/4"°^'")(E ^ 0) (where 4"°^"^) = (m - 1)!/C"\ Eq.(l60]), are the moments 
away from the anomaly): 

Rm = trr'^u r^'^ cos'-{<i>)u-~' (^,0)]^ (135) 

r(m)r(m — 1) z/(0) JO Jo 

Equation Eq. (11351) with i^o{0)/u{0) taken from Eq. (11341) and (u,0) given by 
Eq. fll30p is parameter- free. 

First of all we check that Rm.~>i = 1 as normalization of wave functions requires. 
Applying to Eq. (ll35p the relation (with 6 = m — 1): 

roc 

/ dx x-^+^ f{x) = r ^ /(O) + 0(1) ,5 = m-1^0, (136) 
Jo 

and using Eq. (11341) one immediately obtains /^"^ = = 1. 

The moments Im with m > 0, m 7^ 1 are essentially governed by the w-dependence 
of the generating function $^""^(m, 0). We evaluated numerically the reduced moments 
Rm up to m = 10. The results are given in Fig. 2 and Table 1. 

One can see that the behavior of reduced moments Rm with relatively small m 
{2 < m < 6), Rm ^ i^o/Cn)"'~\ is, indeed, compatible with Eq.([l32]). The best 
exponential fit gives the same tendency of increasing the localization length at the 
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1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1.0000 


0.8347 


0.6703 


0.5321 


0.4252 


0.3467 


0.2908 


0.2519 


0.2255 


0.2083 



Table 1. Reduced moments (m = 1, 2. ..10) at the £' = anomaly for an infinite 
chain in the limit of weak disorder. 



Figure 4. (color line)Reduced moments Rm (red points) in the log-linear scale, 
dashed hne is the exponential fit i?,„ ^ (^o/^an)""^ with CJ/fo = 1-252. 



The 



anomaly as in Eq. fll3ip . However, at the £" = anomaly the extrinsic and intrinsic 
localization lengths are no longer equal to each other, as it is the case at energies away 
from the anomaly: 

C ^ 1.2524, C ^ l-094£o. (137) 

Even more interesting phenomenon takes place for large moments. At m > 6 one can see 
a significant enhancement of the moments compared to their value extrapolated from 
the exponential dependence of at small m. A possible physical meaning of these 
new regime is discussed in a short publication [29] . 



8. Conclusions and open problems 



Eq.f ll30p is the main result of the paper. It gives an exact and unique stationary (i.e. 
site- independent) solution (in quadratures) to the partial differential equation Eq.l lHTj) 
for the generating function $(m, 0) at the center of the energy band {k = 7r/2, E{k) = 0) 
of a weakly disordered chain. The variables u and (j) are associated with slowly varying 
(squared) amphtude and phase of wave functions. The generating function we obtained 
can be used to compute all local statistics of normalized eigenfunctions in the one- 
dimensional Anderson model in the bulk of a long (A^ ^ i) chain. The solution of 
this problem goes beyond the known problem of the Lyapunov exponent and related 
quantities (e.g. density of states and conductance) [HI HSl [171 [IE] , which are completely 
determined by the distribution function of phase V{4>). As we have shown, V{(f)) is 
a descender of the generating function $(m,0) and is related with it in a simple way: 
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P(0) = = 0,0) (see Eq. llSU]) ). Another important application of Vi^cj)) is that 
P(0 = 9/2)/2 is the distribution function Prcf(^) of a random phase 9 of the coefficient 
of reflection from a semi-infinite disordered chain. 

The integrability of the partial differential equation Eq. (l82l) for the generating 
function 0) which we discovered, is a remarkable evidence of a hidden symmetry 
of the problem sX k = 7r/2. Although in the course of derivation we mentioned about 
the sl/2 algebra of operators Eq. flHTl) . which are the building blocks for Eq. fl86l) . and 
about even more extended algebra Eq.f lQOI) . we did not exploit this algebraic content 
explicitly. What would be highly useful is to find the symmetry transformation which 
enables the discovered separation of variables. Moreover, the hidden symmetry which 
survives violation of the chiral (or sublattice) symmetry [281 12] by on-site disorder, could 
be important for other systems (like edge states in Quantum Hall effect or in topological 
insulators) where the chiral symmetry is broken by disorder. Speculating on its nature 
we may surmise that this symmetry might be more naturally formulated in the three 
dimensional space rather than in the two-dimensional space (^, rj) and that it may have 
something to do with the symmetry of the 3d harmonic oscillator. This conjecture is 
fed by an analogy between our main result Eq. f ll30p and the expression for the Green's 
function of the 3d harmonic oscillator problem [26] . The analogy concerns the parameter 
(A in our problem and k in Ref.[26]) entering both the argument and the first index of 
the Whittaker functions in a mutually reciprocal way, as well as the second index of 
the Whittaker functions being ^ in both cases (showing a connection with the parabolic 
cylinder function). Establishing this symmetry would also be useful for studying higher 
order anomalies at /c = nm/n with n > 2 (for the preliminary analysis of the "devil's 
staircase" of these higher-order anomalies, see Ref. [27]). 

We would like to note that the studied anomalies are inherent not only to the 
problem of a disordered chain but might occur in other physical situations where there 
is a periodic perturbation of random amplitude. We mention an analogy between the 
Id localization and the classical system of kicked oscillator studied recently in Ref.|16]. 
According to this analogy the energy-dependent de-Broglie wavelength A^; of a particle 
on a chain is encoded in the frequency of the oscillator while the lattice constant a 
determines the period of the (5-function pulses of the external force ("kicks"), their 
amplitude being proportional to disorder. From this point of view, statistical anomalies 
arise due to sequences of several kicks with amplitudes correlated in time. Correlated 
amplitudes of kicks correspond to exclusive configurations of the local disorder, hence 
the anomalies are weak (for weak disorder) and exhibit themselves only in narrow 
windows around selected energies. Remarkably, the variables ^ and r/, which allowed us 
to factorize the equation for the generating function of the Anderson model, play a role 
of the co-ordinate and the momentum of the kicked oscillator. 

Finally, we applied the exact solution for the generating function to analyze the 
behavior of moments of the eigenfunction distribution at the center-of-band anomaly. 
We have found that relatively small moments behave similar to those outside the 
anomalous region but with a renormalized localization length — >■ t^^, while the 
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larger moments deviate significantly from this dependence. This fact together with 
the appreciable enhancement of the "intrinsic" length ^ 1.252 £o with respect to 
the "extrinsic" length (inverse Lyapunov exponent) ~ 1.094 £o, implies a significant 
change of the form of the "average" eigenfunction at the center-of-band anomaly and 
simultaneously a failure of one-parameter description of eigenfunction statistics. 
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Appendix A. Normalization of $('^")(u = 0,</)): explicit check of the relation 

${'^")(n = 0,0) = P("")(0) for Eqs. (fT30l) and ( l84l) . 



In order to take the limit n — j- in Eq.f ll30p . we re-scale the integration variable 
A — mA and use the identity Wo^u{z) = ^J^K^(^). Introducing a new integration 
variable x = 1/ (8A), we arrive at the following expression for $("")(« = 0, 0): 



4 



$(-)(n = 0, 0) = -^^v'l sin (20)1 (A.l) 



where (see |24 



POO 

/ y/xdxK 1 {ex cos^ (j))Ki {ex sin"^ (f)) (A. 2) 

Jo * 

\/27r^z sin 2 6 3 3 , . / . „x 



COSi 



The latter expression is rather complicated. To find its real part one has to use nontrivial 
identities for the hypergeometric function. Obviously this brute force approach does not 
exploit efficiently the symmetry of the problem. 

It is more advantageous and instructive to exploit the symmetry and perform the 
integration in Eq.( lA.2l) in several elementary steps. Using the representation ( |24] ) 



3' ' V2 



, 1/4 r I - 

K,{z) = (-) I dt e~''°'^Wsmht, (A.4) 



r (!) 



one can perform an elementary integration over x in Eq.( lA.2l) arriving at 



3?/ = — 7=^1 sin COS / (iti (it2vsinhti sinht2 
2v27r Jo ^ 



[e cos^ cosh ti + esin^ cosh ^2]^ ^ ^ 
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Using the identity 

^= r 2^"'^ 212 (A-6) 

[ex + eyj^ [x^ + y^J^ 

(for real a; and y) and introducing new variables: 

yi = cos^ sinhti ; y2 = sin^ sinht2 , 

we obtain: 



m = r ^^v/^/f^^^^ . (A.8) 

TT J I sin (20)1 ^0 [ cos4 + sm^ + 2/f + ^aT 



Re-scaling the variables ?/i(2) = y^g) V *^ ~'~ ^^^^ *^ introducing the polar 

coordinates p and a G (0,7r/2) as y[ = pcosa and 1/2 = psina, we get 



m= (A.9) 

TT J I sin (20) I J cos^ + sin^ 



Ia= Vsina cos a = -E(-, -) = ; (A.IO) 



where 



r = / ^ ^ = -. A.ll 

Jo [l + p2]2 4 ^ ^ 

Collecting things together we arrive at the following final expression for $(«")(n = 0,0), 
Eq.dAl]),: 

$(-)(n = 0, 0) = ^ ^ = , ^ . (A.12) 

n(i JcosV+sinV r2 i ^3 + cos (40) 



This expression coincides with the anomalous probability distribution of phase (0) 
Eq.l lMl) and thus proves the correct choice of the numerical pre-factor in Eq.f llSOp . 
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